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Air  Campaign  and  High-Energy  Laser 
Propagation  Analyses 


Introduction 


This  report  covers  work  LMl  did  from  July  2003  through  September  2004  for  the 
Tactical  Air  Forces  Division  of  Program  Analysis  and  Evaluation  (PA«S:E 
TACAIR)  in  the  OtTice  of  the  Secretary  of  Defense.  The  purpose  for  this  study 
and  its  predecessor  studies  was  to  develop  tools  for  analyzing  tactical  air  forces 
that  meet  three  criteria; 

♦  The  time  required  to  set  the  tools  up  and  run  them  to  generate  usable  in¬ 
formation  should  be  relatively  short — minutes,  not  days.  This  criterion  is 
important  because  analysts  in  PA&E  TACAIR  must  deal  with  a  great 
many  hurried  requests  for  substantial  amounts  of  accurate  information. 

♦  The  tools  should  operate  on  widely  available  personal  computer  (PC)  plat¬ 
forms  without  highly  specialized  software.  This  criterion  is  important  be¬ 
cause  the  analysts  should  be  able  to  use  the  tools  in  readily  available 
environments  without  extensive  training. 

♦  The  tools  should  account  rationally  and  quantitatively  for  the  uncertainties 
inherent  in  air  combat.  This  criterion  is  important  because  PA&E  TAC 
AIR  advises  the  senior  leadership  of  the  Department  of  Defense,  and  these 
clients  need  quantitative  measures  of  both  central  tendencies  and  disper¬ 
sions  to  make  good  decisions. 

To  meet  the  three  eriteria,  LMl  researchers  invested  considerable  effort  in  identi¬ 
fying  the  most  salient  features  of  the  aspects  of  tactical  air  engagements  that  we 
treated  and  in  constructing  models  that  deal  with  them  as  effectively  as  possible 
under  the  limitations  imposed  by  the  first  two  criteria. 

The  report  has  three  major  sections.  The  first,  “An  Application  of  a  Stochastic 
Model  of  an  Air  Superiority  Campaign.”  discusses  details  of,  and  results  from,  the 
Stochastic  Lanchester  Air-to-Air  Campaign  Model  (SLAACM).  This  section  con¬ 
tinues  work  previously  reported  by  LMl  [1  and  2]. 

The  second  section  of  the  report,  ’’Analysis  of  Force  Concentrations  in  Determi¬ 
nistic  Lanchester  Campaigns,”  documents  an  analysis  of  conditions  under  which  a 
weaker  force  can  divide  and  overcome  a  stronger  force. 

The  third  section  of  the  report,  “Analyzing  Laser  Weapons  for  Tactical  Aircraft,” 
discusses  the  algorithms  and  data  being  used  to  analyze  the  potential  performance  of 
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tactical  high  energy  lasers.  It  also  discusses  the  progress  of  transforming  these  analy¬ 
ses  into  a  spreadsheet-based  Airborne  Laser  Infrared  Transmission  (ALIRT)  model. 

Other  support  efforts  provided  under  this  task,  including  a  user's  guide  for  the 
current  version  of  SLAACM  and  a  feasibility  study  for  collection  and  analysis  of 
Battle-of-Britain  air  engagement  and  campaign  data  are  documented  separately. 
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An  Application  of  a  Stochastic  Model 
OF  AN  Air  Superiority  Campaign 


This  section  describes  an  update  of  the  campaign  model  described  in  Stochastic  Models 
of  Air  Superiority  Engagements  and  Campaigns  [1  ]  and  gives  an  example  application. 

The  update  accommodates  impacts  of  varying  air-to-ground  weapon  payload  quantity  and 
quality  of  the  bomber  types  and  the  effect  of  the  defenders’  ability  to  identify  the  aircraft 
types  making  up  an  attack  package.  The  application  treats  a  campaign  scenario  in  which 
an  attacking  force  of  many,  relatively  weak  aircraft  confronts  a  defending  force  of  fewer, 
much  stronger  machines.  In  previous  work,  we  analyzed  individual  engagements  and  the 
combination  of  engagements  into  a  campaign.  Below,  we  address  evaluation  of  the  air-to- 
ground  campaign,  considering  optimal  operations  of  both  attacker  and  defender. 

Scoring  Air-to-G round  Success 

We  wish  to  score  the  value  of  an  attack  package  to  Red  with  a  measure  of  the  value  of  the 
targets  that  the  package’s  payload  is  expected  to  destroy.  First,  let  us  consider  how  the 
value  of  a  target  may  be  expected  to  scale  with  the  yield  of  a  bomb  required  to  destroy  it. 

Value  and  Lethal  Radius  as  Functions  of  Yield 

To  see  how  value  should  vary  as  a  function  of  the  explosive  yield  required  to  de¬ 
stroy  it.  we  considered  the  cost  of  hardening  a  target  against  a  given  yield.  We 
then  assume  that  the  target's  value  is  proportional  to  that  cost. 

To  see  how  the  yield  required  to  destroy  a  structure  varies  with  properties  of  the 
structure,  we  considered  a  simple  model  reinforcing  structure,  shown  in  Figure  1 . 

The  structure  has  depth  (into  the  page)  of  w. 

Figure  1.  Mode!  Structure 

-A 


h 


d 

Pressure  p  on  the  structure’s  leftmost  face  will  induce  both  axial  and  shear  loads  in  the 
reinforcing  columns.  The  axial  load  will  be  pwh“/d,  compressive  on  the  rightmost  col¬ 
umn  and  tensile  on  the  other.  The  shear  loads  on  each  column  will  be  pwh/2.  This  leads 
to  additional  stress  of  maximum  amplitude  close  to  ph7dA.  If  this  stress  is  reacted  by 
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material  of  maximum  stress  Omx,  then  the  thickness  A  of  the  reinforcing  columns  must 
be  ph7d  Omx-  We  assume  that  the  cost  of  reinforcement  is  proportional  to  the  columns' 
volume.  Thus  we  are  led  to  the  conclusion  that  the  cost  of  hardening  the  structure 
against  overpressure  p  is  proportional  to  pwhV  d  Omx,  for  a  given  structure  and  reinforc¬ 
ing  material — that  is,  for  given  w,  h.  d.  and  Omx —  the  cost  of  hardening  is  simply  pro¬ 
portional  to  overpressure  p. 

We  used  dimensional  analysis  to  see  how  pressure  load  varies  with  explosive  yield 
E,  assuming  that  the  pressure  is  dynamic  pressure  induced  by  the  air  velocity  v  im¬ 
mediately  behind  the  blast’s  shock  wave.  In  addition  to  E,  parameters  for  the  blast 
wave  are  initial  air  density  po  and  the  time.  It  follows  that  v  is  proportional  to 
(E/  Since  the  radius  of  the  blast  wave  increases  monotonically  with  time, 

evaluating  v  at  a  given  radius  is  equivalent  to  evaluating  v  at  a  given  time.  Thus  the 
value  of  V  at  a  given  radius  r  is  proportional  to  E*'^'^,  and  the  dynamic  pressure  p  oc 
V*  oc  E  .  Thus  we  are  led  to  take  the  cost  of  hardening  a  structure  as  proportional  to 
the  2/5  power  of  the  explosive  yield  against  which  the  structure  is  hardened. 

Dimensional  analysis  also  lets  us  describe  the  variation  of  the  lethal  radius  n  of 
an  explosion  of  yield  E  against  a  target  hardened  to  withstand  dynamic  pressure  p. 

If  n  oc  E®  p^  then  a  =  1/3  and  b  =  -1/3.  That  is,  rt  (E/p)''^^. 

Value  of  Targets  Destroyed  by  a  Given  Payload 

The  single  shot  kill  probability  (sspk)  that  a  target  is  destroyed  in  one  attack  by  an 
explosion  of  lethal  radius  n  when  delivered  by  a  system  with  a  specified  circular 
error  probable  (CEP)  is 


sspk  =  1  - 


,  yr/cEP- 

2. 


[Eq. 1] 


Recalling  the  variation  of  rt  with  E  and  p,  we  see  that  if  we  refer  tl  to  some  stan¬ 
dard  value  rtrets  hardness  p  to  some  standard  value  pa-r,  and  CEP  to  some  standard 
value  CEPref,  then  ( 1 )  and  the  fact  that  n  oc  (E/p)''^^  imply  that 


sspk  =  1  -  — 


.  ’’Lref  1 

'cep„,Y 

CEP',  1 

P  J  ' 

X  CEP  J 

[Eq. 2] 

Equation  2  gives  us  means  of  determining  sspk  for  varying  yield,  hardness,  and 
CEP.  We  will  assume  that  a  given  aircraft  and  load  combination  will  be  sent 
against  targets  such  that  riat  =  CEPref.  We  further  assume  that  the  yield  is 
matched  to  the  target,  so  that  E/Eref  =  Pret^'p  =  I  •  Then 


sspk  =  1  - 


(-] 

.2v 


CEP^ 

CEP 


[Eq.31 
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If  an  aircraft  with  payload  P  carries  bombs  of  weight  500w,  where  w  is  I,  2,  or  4, 
then  the  number  of  bombs  carried  is 


N  = 


P 

500w 


[Eq.  4] 


Note  that  the  brackets  that  close  only  on  the  bottom  indicate  that  the  result  is  rounded 
down  to  the  next  lowest  integer.  This  is  called  a  “floor  function”  and  is  equivalent  to 
a  truncation  function.  Conversely,  brackets  closed  only  at  the  top  designate  a  “ceiling 
function”  where  the  result  is  rounded  up  to  the  next  higher  integer. 

If  each  bomb  carried  is  dropped  on  a  different  target,  the  expected  number  of  kills  is 


<  targets  killed  >= 


P 

500w 


1- 


CEI^ 
I  ^  CEP- 


[Eq. 5] 


The  value  of  the  targets  varies  as  the  2/5  power  of  the  explosive  yield  required  to 
destroy  them.  In  view  of  this,  the  expected  value  of  the  targets  killed,  relative  to 
some  reference  value,  valuerer,  is 


< 


value 

value 


>= 


P 

500w 


ri 


CEP- 


ref 


CEP- 


[Eq.  6] 


We  used  Equation  6  to  determine  the  values  of  bombers  having  various  values  of 
payload  and  CEP.  Table  I  shows  the  results.  For  completeness,  we  included  fight¬ 
ers,  which  could  be  used  as  bombers. 

Table  1.  Aircraft  Values  as  Bombers 


Aircraft 

Weapon  parameters 

Weapon  value  versus: 

Summary  results 

Payload 

CEP 

SSPK 

500-lb 

target 

1,000-lb 

target 

2,000-lb 

target 

Average 

value 

Relative 

average 

value 

RED F1 

1,000 

1 

0.50 

1.0 

0.7 

0.0 

0.6 

1.0 

RED F2 

1,000 

1 

0.50 

1.0 

0.7 

0.0 

0.6 

1.0 

RED F3 

2,000 

1 

0.50 

2.0 

1.3 

0.9 

1.4 

3.0 

RED F4 

3,000 

0.6 

0.85 

5.1 

3.4 

2.2 

3.6 

6.0 

RED F5 

6,000 

0.1 

1.00 

12.0 

7.9 

5.2 

8.4 

15.0 

RED FB1 

2,000 

0.6 

0.85 

3.4 

2.3 

1.5 

2.4 

4.0 

RED B2 

6,000 

0.4 

0.99 

11.8 

7.8 

5.2 

8.3 

15.0 

RED_B3 

18,000 

0.1 

1.00 

36.0 

23.8 

15.7 

25.1 

45.0 
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Scoring  of  Aircraft 

Value  to  Red  of  a  Blue  Fighter  Killed 

Assuming  Blue  aircraft  have  8  missiles,  with  sspk  of  0.8,  sending  a  single  flight  of  4 
Blues  against  two,  12-aircraft  attack  packages  gives  less  than  83  percent  probability 
of  killing  all  24  opponents.  Therefore,  Blue  is  likely  to  use  its  aircraft  in  flights  of  4 
against  only  one  attack  package  of  12  aircraft.  In  one  mission,  each  Blue  aircraft  is, 
therefore,  “worth”  3  Red  aircraft.  The  average  value  of  Red  aircraft  is  roughly  6,  and 
so  we  took  the  value  to  Red  of  killing  a  Blue  aircraft  to  be  1 8. 

With  the  scoring  values  thus  developed,  we  took  the  Red  payoff  of  an  attack 
package  to  be  four  times  the  product  of  the  probability  that  the  bombers  got  past 
Blue  defenders  and  the  value  of  the  individual  bombers  as  shown  in  Table  1,  plus 
1 8  times  the  expected  number  of  Blue  aircraft  killed. 

We  developed  distributions  of  Red  and  Blue  losses,  and  the  probability  that  the 
bombers  got  through  the  defenders  (that  is,  the  probability  that  not  all  8  escorts 
were  killed),  from  the  following  variation  of  the  classic  stochastic  Lanchester  en¬ 
gagement  model:  first,  the  4  advanced  escorts  engaged  the  Blue  defenders  in  a 
4-Vs-4  engagement  to  annihilation  on  both  sides.  Then,  any  surviving  Blue  air¬ 
craft  engaged  the  4  close-in  escorts;  again  all  engagements  were  to  annihilation.' 

Values  for  Blue 

We  assigned  Blue’s  penalty  for  losing  a  Blue  fighter  as  40.  We  assigned  Blue's 
benefit  for  killing  a  Red  bomber  as  that  bomber’s  value  to  Red. 

Campaigns  and  Their  Results 

We  modeled  campaigns  in  the  following  way.  Starting  with  the  Red  order  of  bat¬ 
tle,  we  developed  an  optimal  Red  assault  as  the  set  of  assault  packages  that 
maximized  Red’s  expected  payoff,  subject  to  the  constraints  imposed  by  the  num¬ 
bers  of  aircraft  of  the  various  types  available.  An  assault  package  is  described  by 
the  ordered  triple  [advanced  escort,  close  escort,  bomber].  There  are  45  possible 
package  options,  if  Red  uses  only  Red  FBl,  Red_B2,  and  Red_B3  as  bombers. 

We  modeled  Blue’s  response  in  two  ways:  in  the  first,  we  assumed  that  Blue  had  perfect 
knowledge  of  the  aircraft  types  making  up  each  attack  package.  I  hat  enabled  Blue  to  re¬ 
spond  by  attacking  that  set  of  Red  attack  packages  that  maximized  Blue’s  expected  payoff. 

In  the  second,  we  assumed  that  Blue  had  no  information  about  the  aircraft  in  an 
attack  package  before  attacking  it.  In  this  case,  the  Red  attack  packages  are  inter¬ 
cepted  entirely  at  random. 


'  “Annihilation” 
destroyed. 


here  means  that  the  combat  runs  until  all  the  aircraft  of  one  of  the  sides  are 
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The  following  sections  show  day-by-day  results  for  the  two  cases  of  Blue’s 
knowledge  of  the  composition  of  Red’s  attack  packages. 


Results  When  Blue  Has  Adequate  Battle  Management  and  Perfect 
Knowledge  of  Red’s  Packages 


We  first  treat  a  campaign  in  which  Blue  has  perfect  knowledge  of  the  aircraft 
types  in  each  attack  package  and  has  adequate  battle  management  capability  to 
assign  the  M  available  Blue  flights  to  any  M  of  the  attack  packages. 


On  the  first  day.  Red’s  optimization  led  to  the  set  of  attack  packages  shown  in  Table  2. 
Table  2.  Day  One  Red  Attack  Packages 


Number  of 

Advanced 

Close 

Bomber 

packages 

escort 

escort 

29 

Red  F5 

Red  F5 

Red B3 

2 

Red  F5 

Red  F5 

Red  B2 

6 

Red  F4 

Red  F4 

Red  FBI 

28 

Red  F3 

Red  F3 

Red  FBI 

50 

Red  F2 

Red  F2 

Red  FBI 

41 

Red  F2 

Red  F2 

Red  B2 

40 

Red  FI 

Red F1 

Red  FI 

1 

Red  FI 

Red_F1 

Red  B2 

Blue’s  optimal  response  was  to  send  all  25  available  fighter  flights  against  the 
most  threatening  attack  packages,  i.e.,  those  with  Red_F5  for  both  escorts  and 
Red_B3  as  the  bomber. 

We  continued  this  campaign,  taking  e.xpected  values  of  the  losses  of  all  aircraft. 
Figures  2-5  show  how  it  evolved. 

Figure  2.  Bomber  Value  Delivered 
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By  “Bomber  Value  Delivered,”  we  mean  the  total  expected  destructive  value  of  the  bomb¬ 
ers  not  killed  by  Blue  defenders,  either  because  their  flights  were  not  intercepted  (the  most 
common  cause)  or  because  their  package’s  escorts  destroyed  all  four  Blue  defenders. 
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Figure  3.  Red  Fighters  Available 
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Figure  4.  Red  Bombers  Available 
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Figure  5.  Blue  Fighters  Available 
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Results  When  Blue  Encounters  Red  Packages  at  Random 


For  random  interceptions,  if  there  are  J  distinct  types  of  attack  packages,  with  N,  pack- 

J 

ages  of  type  j,  and  a  total  of  N  =  attack  packages  launched  in  the  Red  assault. 

i=i 

and  M  defending  flights  respond,  then  the  probability  P(ni,  n:,  ...,  nj)  that  ni  packages 
of  type  I,  n2  packages  of  type  2,  ...,  nj  packages  of  type  J  are  intercepted  is 


P(n|,n2,-../nj)  = 


M!  (N-M)!  J  Nj! 

On.! 


[Eq.  7] 


The  distribution  of  the  attack  packages  intercepted  when  Blue  intercepts  them  entirely 
at  random  is  dispersed,  in  the  sense  that  many  distinct  sets  of  intercepted  packages, 
characterized  by  the  n,  of  Equation  7,  have  nearly  the  same  probability.  The  impact  of 
this  dispersion  may  be  mitigated,  however,  by  the  fact  that  the  sets  of  packages  that  are 
intercepted  with  relatively  high  probability  all  have  nearly  the  same  bomber-damage 
value  and  blue  losses.  Figure  6  illustrates  this,  showing  probability,  expected  bomber- 
damage  value,  and  expected  total  Blue  losses  for  the  12  most  likely  interceptions. 

Figure  6.  Interception  Statistics 


□  Probability  x  10,000 
■  Damage  x  0.001 

□  <Total  Blue  Loss> 


To  determine  damage,  and  the  next  day's  Blue  force,  we  assumed  that  the  mode 
of  the  Equation  7  distribution  occurred.  In  all  the  cases  that  we  considered,  the 
mode  did  not  differ  significantly  from  the  mean,  rounded  to  integer  values. 

Figures  7  and  8  compare  damage  value  and  Blue  aircraft  available  for  the  first  two 
days  of  the  "“perfect  information”  and  ""random  intercept”  campaigns.  Figure  7 
shows  the  substantial  value  of  Blue  being  able  to  put  its  assets  against  only  the 
most  threatening  attack  packages.  The  total  damage  values  are  34,100  units  for 
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the  ■"random”  case  and  12,152  damage  units  for  the  "perfect”  case.  Perfect  knowl¬ 
edge,  thus,  results  in  a  64  percent  reduction  in  enemy  damage.* 

Figure  7.  Bomber  Damage  Value  Comparison 
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Perfect 


— Random 


Figure  8  shows  that  Blue’s  initial  losses  will  be  less  in  the  "random”  case,  because 
Blue  does  not  fight  as  many  engagements  with  Red’s  best  fighters  as  in  the  "perfect” 
case.  We  see  here  that,  in  order  to  achieve  the  dramatic  reduction  in  Red  damage 
with  perfect  knowledge.  Blue  must  not  flinch  in  the  face  of  severe  initial  losses. 

Figure  8.  Blue  Aircraft  Available  Comparison 


^  The  corresponding  values  in  tons  of  bombs  dropped  are  7,076  for  the  random  case  and  2,768 
for  the  perfect  case,  resulting  in  a  61  percent  reduction  due  to  perfect  knowledge. 
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Figure  9  shows  the  variation  of  the  numbers  of  the  various  types  of  Red  aircraft. 

Figure  9.  Red  Aircraft  Available 
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Figure  9  shows  consistent,  linear  decreases  in  the  numbers  of  each  type  of  Red  air¬ 
craft.  This  contrasts  sharply  with  the  variations  shown  in  Figures  3  and  4.  For  the 
campaign  described  in  those  figures.  Blue’s  perfect  information  about  the  composi¬ 
tion  of  Red  packages  allows  Blue  to  concentrate  systematically  on  the  most  threat¬ 
ening  formations.  With  the  random  interceptions,  and  the  great  dominance  of  Blue’s 
aircraft,  the  rate  of  decrease  of  a  Red  aircraft  type  will  be  roughly  proportional  to 
the  fraction  of  the  total  number  of  Reds  that  the  type  represents.  Those  fractions  are 
roughly  constant  through  the  campaign’s  first  4  days. 


Analysis  of  Force  Concentration  in 
Deterministic  Lanchester  Campaigns 


In  this  section  we  show  that,  in  the  context  of  the  force-on-force  Lanchester 
model,  the  smallest  number  of  engagements  with  which  an  inferior  force  can  an- 

I 


nihilate  a  stronger  one  by  repeated  force  concentration  is 


Po 


,  where  po  is  the 


initial  force  ratio.  We  also  derive  some  optimal  force-concentration  strategies. 


Background 

The  analysis  of  force  concentration  is  a  major  theme  of  Lanchester's  seminal  pa¬ 
per  [3].  Subsequent  work  has  extended  Lanchester’s  results,  for  example,  by  con¬ 
sidering  optimal  allocations  of  elements  of  a  superior  force  [4]. 


In  the  context  of  the  Lanchester  force-on-force  model,  which  is  defined  in  Taylor 
[5],  a  weaker  force  can  defeat  a  stronger  one  by  force  concentration  in  a  series  of 
engagements.  Lanchester  [3]  cites  Nelson's  plan  for  the  battle  of  Trafalgar  as  an 
e.xample.  The  British  fleet  was  weaker  than  the  French-Spanish  combined  fleet  that 
they  faced,  by  a  force  ratio  of  about  3:4.  In  his  plan.  Nelson  concentrated  essentially 
all  his  force  on  half  the  enemy  forces.  The  British  would  have  had  a  roughly  2-1 
force-ratio  advantage  in  this  initial  engagement,  which  they  would  have  won  with 
sufficient  survivors  to  defeat  the  other  half  of  the  enemy  force.  Thus  Nelson's  plan 
would  have  been  expected  to  result  in  victory  in  a  sequence  of  two  engagements. 

This  suggests  the  possibility  of  a  general  strategy  for  a  weaker  force  to  annihilate  a 
stronger  one,  by  successively  concentrating  all  available  force  on  selected  fractions 
of  the  enemy  force.  Let  us  call  the  weaker  force  Red,  the  stronger.  Blue.  Perhaps  if 
Red  concentrates  all  its  force  on  a  fraction  fo  of  the  Blue  force.  Red  can  win  that 
engagement  and  also  a  subsequent  engagement  between  the  survivors.  If  so.  Red 
will  win  the  campaign  in  two  engagements.  If  that  is  not  possible,  perhaps  by  con¬ 
centrating  all  its  forces  on  a  fraction  fo  of  the  Blue  force,  and  then  concentrating  the 
surviving  Red  force  on  a  fraction  f|  of  the  surviving  Blue  force.  Red  can,  in  a  third 
engagement,  annihilate  all  the  Blues  who  survive  the  second  engagement,  and  so 
win  the  campaign  in  three  engagements.  If  that  is  not  possible,  perhaps  in  a  similar 
way  Red  can  win  in  four  engagements,  or,  if  not  in  four,  perhaps  in  five,  and  so  on. 

The  purpose  of  this  analysis  is  to  examine  this  strategy  in  general.  We  will  exhibit  a 
lower  bound  on  the  number  of  such  engagements  and  derive  some  optimal  strategies 
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statement  of  the  Problem 


A  Red  force  with  Ro  units  of  effectiveness  kr  confronts  a  Blue  force  with  Bo  ele¬ 
ments  of  effectiveness  kb.  The  Red  force  is  weaker,  in  the  sense  that  force  ratio  po 
of  these  initial  forces. 


[Eq. 8] 


is  less  than  one.^  Thus  Red  will  lose  an  engagement  with  the  entire  Blue  force. 


But  perhaps  if  Red  concentrates  all  his  force  on  a  fraction  fo  of  the  Blue  force,  he 
can  win  that  engagement  and  also  a  subsequent  engagement  between  the  Red  and 
Blue  forces  that  survive  it.  If  so.  Red  will  win  the  campaign  in  two  engagements. 
If  that  is  not  possible,  perhaps  by  concentrating  all  his  forces  on  a  fraction  fo  of 
the  Blue  force,  and  then  concentrating  his  surviving  force  on  a  fraction  f|  of  the 
surviving  Blue  force.  Red  can,  in  a  third  engagement,  annihilate  all  the  Blues  who 
survive  the  second  engagement,  and  so  win  the  campaign  in  three  engagements. 

The  problem  we  consider  is  to  determine  the  smallest  number  of  engagements 
with  which  Red  can,  by  concentrating  its  force  on  specified  fractions  fo,  f|,  ...  of 
the  Blue  force  confronting  him,  annihilate  the  Blue  force.  We  also  wish  to  deter¬ 
mine  choices  of  sequences  fo,  fi,  ...  that  minimize  Red’s  losses. 


Results 


Our  principal  results  can  be  stated  conveniently  as  two  theorems: 

1 .  Theorem  I .  A  necessary  and  sufficient  condition  that  Red  ean  win  the 
campaign  in  as  few  as  N  engagements,  N  =  2,  3,  ....  is  that  po  >  l/N. 

2.  Theorem  2.  If  Red  wins  the  campaign  in  N  engagements,  the  sequence  fo, 
fi,  ....  fN.2  that  minimizes  Red’s  losses  is 


[Eq.  9] 


and  ,  the  square  of  the  number  of  surviving  Reds,  is  given  by 


[Eq. 10] 


^  kxX"  is  also  referred  to  as  the  fighting  strength  of  the  forces,  and  po  referred  to  as  the 
strength  ratio. 
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where  =  l/N . 

Two  lemmas  facilitate  proofs  of  the  theorems: 

1 .  Lemma  1 .  The  function  gN(x), 

gN(x)  =  -^('-x)^ [Eq.  Ill 
has  a  unique  minimum  at  x  =  1/(N  +  1 ). 

Proof:  g'tsj  (x)  is  continuous,  negative  for  x  <  1  /(N+ 1 )  and  positive  for  x  >  I  /(N+ 1 ). 

It  will  also  be  useful  to  note  that  gN(l/(N  +  I ) )  =  l/(N  +  1 ). 

2.  Lemma  2.  A  necessary  and  sufficient  condition  that  there  are  values  of  f  in 
(0, 1 )  such  that 


p-f- 


[Eq. 12] 


is  that  p  >  1/(N  +  I ),  for  integer  N  >  2. 

Proof:  Inequality  12  implies 

p>:l(l-f)-+f‘=gN(0  [Eq.  13] 

Necessity  then  follows  from  Lemma  I .  By  the  continuity  of  gN(x),  if  p  > 
1/(N+1 )  there  is  a  neighborhood  of  f  =  1/(N+I )  for  which  Equation  13, 
and  consequently  Inequality  12,  is  satisfied.  Since  0  <  1/(N  +  l)<  I  forN 
>  2,  there  is  also  a  sub-interval  of  (0,1 )  within  which  Inequality  12  is  met. 
and  that  establishes  sufficiency. 


Proof  of  Theorem  1 :  We  proceed  by  induction  on  N.  Let  N  =  2.  For  Red  to 
win  the  first  engagement. 


_  Po_  ^  I 


[Eq.  14] 


When  Inequality  14  is  satisfied,  the  second  engagement  will  have  force  ratio  pi 
given  by 


R5-k,B5fo^  _  po-f,? 

O-fo)- 


[Eq. 15] 
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Thus  Red  will  annihilate  the  Blue  force  in  two  engagements  if  and  only  if  (iff)  pi 
>  I  and  Inequality  14  is  met.  But  Inequality  14  is  satisfied  for  all  fo  e  (0,  I )  that 
make  pi  >  I,  since 


Po  ^0 


>I^Po  >(l-fo)‘  +fo 


(l-fp)- 


+  1 


[Eq.  16] 


By  Lemma  2,  there  are  such  values  of  fo  iff  po  >  1/2,  so  the  theorem’s  statement  is 
correct  for  N  =  2. 


Now  suppose  the  statement  is  correct  for  N  =  K.  Then  Red  can  win  in  K  +  I  en¬ 
gagements  if  (6)  is  met,  and  if  it  is  possible  to  choose  fo  such  that  pi  >  l/K.  By  (8) 
and  Lemma  2,  this  is  possible  iff  po  >  1/(K  +  I ),  which  completes  the  induction. 


♦ 


Corollary  to  Theorem  1 ;  The  smallest  number  N  of  engagements  with  which  an 
inferior  force  can  annihilate  a  superior  one  by  force  concentration  in  successive 


engagements  is 


J_ 

Po 


Proof  of  Theorem  2:  Again  we  proceed  by  induction.  First,  we  note  that 
when  Red  annihilates  the  Blue  force  in  a  sequence  of  N  engagements,  the  number 
Rn  of  surviving  Red  elements  is  given  by 


Rl  =  +(l  -/o)V,-  -0-foy-0-A)'f2  +■■ 


+  (I  -  ./o  )■  (1  -  -(I  -  A-3  )-  (I  -  /lV-2  )■  1 


[Eq. 17] 


Thus  the  number  of  Red  survivors  will  be  maximized  by  the  choice  fo,  fi,  ...,  fNo 
that  minimizes  Jn,  where 

./.(/o,/ . /»-,)= U,’+o-/o)'/’+o-/«)'o-y;)V2-+- 

+  (l-/o)'(l -/v-3)'/v-2  [Eq.  18] 

Note  that  Equation  18  implies  a  recursion  of  the  Jk, 

. ,fK-,)  =  fo  +(l-fo)'JK(f.^f2 . fx-i)  [Eq.  19] 
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If  N  =  2,  then  J2(fo)  =  gi(fo)  and.  by  Lemma  I,  the  minimizing  value  of  fo  is  1/2, 
and  Jj,  the  minimum  value  of  Ji,  is  1/2.  Therefore  the  statement  of  Theorem  2  is 
correct  for  N  =  2. 

Now  suppose  the  statement  is  correct  for  N  =  K.  That  is,  suppose  JkCL),  fu  •••,  fk 
.2)  is  minimized  by  fi=  1/(K  -  j),  j  =  0,  1,2,  ...,  K  -  2,  and  thatjJ^,  the  smallest 
value  of  Jk,  is  given  by  =  I  /  K  . 

Referring  to  Equation  19,  we  see  that  the  coefficient  of  Jk(fi.  f:.  •••.  fk-i)  is  posi¬ 
tive  for  all  fo,  so  Jk+i  will  be  minimized  by  choosing  fi,  f2,  ...,  fk-i  to  minimize 
JkCfi.  f:-.  •••■,  fk-iL  and  then  choosing  fo  to  minimize  f^  +  (I  -  f^)' ,  where 

denotes  the  minimum  value  of  Jk-  By  the  hypothesis  that  the  statement  is  true 
for  N  =  K,  we  have  fj  =  l/(K  +  I  -  j),  j  =  1, 2,  ...,  K  -  I .  Also  by  that  hypothesis, 

=  l/K  .  Then,  by  Equation  19,  fo  must  be  chosen  to  minimize  gK(fo)-  By 
Lemma  1,  that  implies  that  fo  =  1/(K  +  I ),  and  that  =  1  /  (K  + 1),  which  com¬ 
pletes  the  induction. 

Discussion 

We  believe  the  results  reveal  an  interesting  feature  of  the  Lanchester  force-on- 
force  model.  We  also  believe  that  the  results  are  useful,  not  necessarily  because 
commanders  may  apply  them  directly,  as  Admiral  Nelson  apparently  did  at  Tra¬ 
falgar,  but  rather,  because  it  will  prove  useful  to  consider  actual  force  concentra¬ 
tion  strategies  in  the  light  of  knowledge  of  the  limits  on  force  concentration,  and 
of  optimal  concentration  strategies. 

Lanchester  [3]  points  out  that  dividing  an  opposing  force  into  two  equal  segments 
gives  the  lowest  possible  total  force.  This  is,  of  course,  an  example  of  Lemma  1 
for  N  =  I . 

While  this  report  was  in  preparation,  we  learned  of  a  paper  by  David  Bitters'*, 
which  obtains  the  principal  result  by  another  method. 


■’  David  Bitters,  “Efficient  Concentration  of  Forces,  or  How  to  Fight  Outnumbered  and  Win." 
Naval  Research  Logistics  42.  pp.  397-418  (1995). 
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Analyzing  Laser  Weapons 
FOR  Tactical  Aircraft 

Background 

There  is  continuing  interest  in  the  potential  for  use  of  airborne,  tactical,  high  en¬ 
ergy  laser  weapons.  Early  efforts  in  this  area  focused  on  10.6  micron  (p)  wave¬ 
length  carbon  dioxide  (CO2)  and  3.8  p  deuterium  fluoride  (DF)  chemical  lasers. 
Current  interest  is  focused  on  1 .3  p  chemical  oxygen  iodine  (COIL),  and  electri¬ 
cally-powered  solid  state  lasers.  The  most  mature  solid  state  laser  technology 
flashlamp-pumped  1 .06  p  neodymium/yttrium  aluminum  garnet  (Nd:YAG),  how¬ 
ever,  diode  pumping,  and  alternative  lasing  elements  and  substrates  are  under  ac¬ 
tive  consideration.^ 

Practical  tactical  laser  weapons  are  conceded  to  be  far  in  the  future,  but  it  is  not 
too  early  to  investigate  potential  laser  performance  and  weapon  utility.  During 
this  year  we  investigated  the  propagation  properties  of  the  iodine  laser  in  the  tac¬ 
tical  environment.  Future  studies  will  include  similar  propagation  studies  of  solid 
state  lasers  and  engineering  investigation  of  both  COIL  and  solid  state  lasers. 

Overview 

In  this  section  we  describe  a  MSExcel  workbook  implementation  of  a  method  for 
calculating  the  peak  fluence  delivered  to  an  airborne  or  ground-based  target  from  a 
laser  weapon  on  an  airborne  platform,  when  the  beam  propagates  through  the  in¬ 
homogeneous,  absorbing,  scattering  atmosphere  (refraction  and  extinction)  which 
is  in  turbulent  motion  (turbulence),  and  when  heating  of  the  air  by  the  beam 
(blooming)  may  have  significant  effects.  Results  reported  here  are  for  one  laser 
wavelength,  that  of  the  Chemical  Oxygen-Iodine  Laser  (COIL),  1 .3 1 52  microns. 
E.xtending  our  methods  to  other,  single-line  wavelengths  is  straightforward. 

The  paragraphs  below  describe  our  methods  for  calculating  the  effects  of  atmos¬ 
pheric  refraction,  extinction,  turbulence,  and  blooming  on  laser  weapons  for  tactical 
aircraft.  Its  following  four  sections  give  details  of  our  treatments  of  these  etTects. 

Extinction 

Our  basic  source  of  extinction  information  was  PcLnWin,  a  graphical  interface  to  the 
FASCODE  atmospheric  code  provided  by  the  Ontar  Corporation.^  PcLnWin  is  fully 
documented  in  such  Ontar  publications  as  PcLnWin  Manual,  Version  J,  October 
1999. 


^  Solid  state  lasers  typically  include  a  lasing  element  such  as  neodymium  dispersed  in  a  host 
structural  material  such  as  yttrium  aluminum  garnet. 

*  Ontar  Corporation,  9  Village  Way,  North  Andover,  Massachusetts  01 845. 
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PcLnWin/FASCODE  requires  certain  user  inputs,  while  providing  default  values 
for  many  FASCODE  inputs.  The  user  inputs  that  we  used  were 

♦  Line  file.  Using  data  provided  with  PcLnWin,  we  built  a  line  file  covering 
the  wavenumber  interval  from  7224.491  cm  '  to  7983.610  cm  ',  roughly 
±5  percent  from  the  COIL  wavenumber  of  7604.523  cm  '.  The  file  in¬ 
cludes  Ni,  O:,  H2O,  CO2,  O3,  N2O,  CO,  and  CH4  molecules. 

♦  Almo.sphere.  We  developed  results  for  three  of  PcLnWin’s  available  at¬ 
mosphere  models:  the  1976  U.  S.  Standard  Atmosphere,  and  two  models 
assembled  by  Ontar  from  various  sources  to  treat  the  winter  atmosphere  in 
midlatitudes,  and  the  summer  atmosphere  in  mid  latitudes.  We  will  refer 
to  these  three  atmosphere  models  as  Std,  MEW,  and  MLS,  respectively. 

♦  Aerosol.  We  chose  FASCODE's  Rural  aerosol.  This  choice  gives  23  km 
visibility  at  the  surface. 

♦  Tropo.sphere/Stratosphere.  This  transition  was  determined  by  the  atmos¬ 
phere  model  used. 

♦  Geometry.  For  our  initial  calculations,  we  used  the  PcLnWin  option  of 
specifying  observer  height,  final  height,  and  earth-centered  angle.  In  sev¬ 
eral  cases  we  verified  that  the  refracted  path  that  we  computed  (our  compu¬ 
tations  are  described  below  in  the  section  “Path  Integrals")  agreed  with  the 
one  computed  by  FASCODE.  Later,  as  described  under  “Blooming,"  we 
developed  analytic  models  of  the  altitude  variations  of  aerosol  absorption 
and  scattering,  and  of  molecular  absorption,  to  treat  STD,  MLW  and  MLS 
atmospheres.  We  also  developed  analytic  models  of  the  altitude  variation  of 
refractive  index  for  the  MLW  and  MLS  atmospheres.  The  STD  atmos¬ 
phere’s  refractive  index  variation  is  analytic,  as  described  below  under 
“Numerical  Evaluation  of  Path  Integrals." 


Turbulence 

Our  principal  parameter  for  turbulence  effects  is  the  Strehl  ratio  SR. 


Maximum  fluence  with  turbulence 
Maximum  fluence  in  diffraction  lim  it 


-f— 1 
71 1  XL  J 


—(Maximum  fluence  with  turbulence) 


[Eq.  20| 


In  Equation  20,  D  is  the  diameter  of  the  laser's  telescope’s  mirror,  and  X  is  the 
wavelength  of  the  radiation.  L  is  the  length  of  the  laser-to-target  path  undisturbed 
by  turbulence,  and  P  is  the  laser’s  power. 

We  used  the  expression  for  SR  given  by  Sasiela  in  Equation  8. 1 3,  on  page  1 64  of  [6]: 
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[Eq.2l] 


SR  = 


24  ^(-1)"  r(3/2  +  5n/6)r(n  +  6/5) 
Syfnh  n!  r(3  +  5n/6)r(n  +  ll/5) 


3.44 


v*'o  y 


We  also  used  the  asymptotic  expansion  for  large  values  of  n/D  given  in  Sasiela's 
Equation  8.15. 

The  parameter  ro  is  the  coherence  diameter  imposed  by  turbulence,  that  is,  the  dis¬ 
tance  perpendicular  to  the  propagation  direction  over  which  turbulent  perturbations  in 
refractive  index  allow  coherence  in  the  electromagnetic  wave  from  the  laser.  We  used 
Sasiela's  expression  for  the  coherence  diameter  of  a  spherical  wave  with  focus  at  the 
target. 


ro,  =  (0.423  [Eq.  22] 

Following  Sasiela’s  Equation  4.6,  we  evaluated  115/3  as  the  path  integral 

P5/3  =  |C‘(z)(L-z)'/Mz  [Eq.23] 

0 

where  denotes  the  refractive  index  structure  parameter  of  the  turbulence. 

Cp  varies  with  height,  geographic  location,  time  of  day,  and  season.  For  cases  in¬ 
volving  propagation  to  the  ground,  we  used  the  Hufnagel-Valley  model  for  C^, 

C;  =8.2xlO'-^W‘h“’e"'’/'  '’ +  2.7xl0’'^e“'’/'-^ +Ae’''/°  '  [Eq.  24] 

The  expression  of  Equation  24  gives  in  m'^^^  for  h  in  km,  W  in  m/sec,  and  A  in 

We  used  W  =  21  m/sec  and  A  =  1.7x10''“’  With  this  choice  of  parameters, 
the  Hufnagel-Valley  model  is  called  the  Hufnagel-Valley  5/7  model,  because  for  ver¬ 
tical  propagation  downward  through  the  entire  atmosphere  it  gives,  in  standard  mod¬ 
els.  a  coherence  radius  of  5  cm  and  an  isoplanatic  angle  of  7  microradians,  for  0.5 
micron  radiation. 

In  some  cases  of  air-to-air  engagement,  we  modeled  as  2xCLEAR  I  Night. 

This  model  is 


log  IQ  =  A  +  Bz  +  Cz"  +  Dexp[{-  0.5(z  -E)/F}“  J 

where  parameters  A,  B,  C,  D,  E,  and  F  are  given  in  Table  3. 
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Table  3.  CLEAR  I  Night  Parameters 


z 

A 

B 

C 

D 

E 

F 

I.23<z<2.l3 

-10.7025 

-4.3507 

0.8141 

0 

NA 

NA 

2.l3<z<l0.34 

-16.2897 

0.0335 

-0.0134 

0 

NA 

NA 

I0.34<z<30 

-17.0577 

-0.0449 

-0.0005 

0.6181 

15.5617 

3.4666 

Numerical  Evaluation  of  Path  Integrals 

We  evaluated  the  path  integral  in  Equation  23,  and  certain  other  path  integrals, 
numerically.  To  do  so,  we  needed  to  compute  the  path  taken  by  the  laser's  output 
through  the  atmosphere  undisturbed  by  turbulence. 

Effects  of  oxygen  and  nitrogen  molecules  dominate  the  atmosphere's  refractive 
index.  Following  the  expression  for  n  given  by  Thomas  and  Duncan  [7],  we  used 

n  =  0.079P/T  [Eq.  25] 

where  P  is  pressure  in  atmospheres  and  T  is  temperature  in  °K. 

For  MLS  and  MEW  atmospheres,  we  used  analytic  models  of  PcLnWin  refractive 
index  outputs.  Calculation  of  the  refracted  ray  requires  derivatives  of  refractive 
index  variation  with  height,  so  a  continuously  differentiable  model  is  desirable. 
We  made  the  models  by  linear  interpolation  in  the  divided  differences  of  refrac¬ 
tive  index  parameter,  and  integrated  the  resulting  piecewise  linear,  continuous 
function  with  an  added,  adjustable  constant  of  integration  to  obtain  a  continuously 
differentiable  model  of  refractive  index. 

Figure  10  shows  the  comparison  between  our  model  of  the  refractive  index  pa¬ 
rameter  le*’(n  -  I ),  where  n  is  refractive  index,  and  the  PcLnWin  data  for  MLS. 

Figure  10.  Verification  of  MLS  Refractive  Index  Model 


--♦—Data 
— Model 
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Figure  10  shows  the  comparison  between  divided  differences  of  the  refractive  in¬ 
dex  parameter,  and  the  derivative  of  our  model  of  the  refractive  index  parameter 
for  MLS.  The  tropopause  is  clearly  noticeable  in  Figure  1 1 . 

Figure  11.  Verification  of  Derivative  of  MLS  Refractive  Index  Model 
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The  refractive  index  variation  of  Thomas  and  Duncan  [6]  is  not  large;  it  is  about 
0.0003  near  the  surface.  Also,  tactical  laser  weapons  will  be  used  at  medium 
ranges,  generally  no  greater  than  100-200  km.  It  is  possible  that  using  straight 
paths  would  not  greatly  distort  integrals  like  that  of  Equation  23. 

Two  considerations  militate  against  neglecting  refraction,  however.  First,  refrac¬ 
tion  is  noticeable  at  extreme  ranges  and  low  weapon  altitudes  ( 

Figure  12).  Second,  and  deciding,  varies  sharply  with  height,  so  that  path 

variations  may  have  significant  effects.  For  e.xample,  differs  by  more  than  10 
percent  on  substantial  parts  of  the  two  paths  shown  in 

Figure  1 2  (please  note  that  the  great  difference  between  the  scales  of  the  ordinate 
and  abscissa  of  that  figure  distorts  geometry). 
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Figure  12.  Refracted  and  Straight  Paths 


I  ■  Earth  - Straight  - Refracted  | 


Accordingly,  we  determined  refracted  paths  numerically.  We  did  so  by  determin¬ 
ing  the  path  of  minimum  propagation  time  from  weapon  to  target.  We  used  the 
coordinate  system  shown  in  Figure  1 3 

Figure  13.  Coordinate  System 


Target 


Specifically,  the  variational  problem  is 


Sl 

min  J 
h(S) 


ylih')-  +{\  +  h/R)-dS 
C(h(S)) 


|Eq.  26] 


where  R  is  the  earth’s  radius  and  C  is  the  speed  of  light.  The  familiar  condition 
for  the  minimizing  h(s). 
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af  d  df 

ah  ds  ah' 


where  f(S,h.h')  is  the  integrand  of  Equation  26,  leads  to  the  system 
^^,_v(S)C(l  +  h(S)/R) 


,_Vl-v-C-  l  +  h/R  I  ac 
RC  Vl-v-C“  C-  ah 


[Eq. 27] 


of  ordinary  differential  equations.  This  system,  with  the  boundary  conditions 

h(0)  =  hw,  h(SL)  =  hr  [Eq.  28] 


determines  h(s). 

A  Consistency  Check  Integral 

Our  model  for  the  Strehl  ratio  rests  on  the  method  of  smooth  perturbations  (MSP)  for 
treating  the  nonlinear  Riccati  equation  for  the  logarithm  of  a  field  component,  that 
arises  in  Rytov's  method  for  determining  the  effects  of  turbulence.  Sasiela  [5]  Section 
2.2,  pp.  30  et  seq.,  notes  that  a  necessary  condition  for  the  consistency  of  this  method  is 

ky‘JC;(h(z))(L-z)*/‘dz<l 

0 

We  evaluated  this  integral  numerically  to  check  on  the  consistency  of  the  MSP. 

Numerical  Method 

We  determined  contours  of  constant  fluence,  accounting  for  atmospheric  extinc¬ 
tion  and  turbulence,  in  this  way:  for  a  given  value  of  desired  fluence  Id,  and  for 
fixed  laser  height  hw  and  target  height  hj,  we  varied  earth-surface  distance  S  to 
find  values  SI  and  S2  that  bracketed  the  value  of  S  that  made 


ld,ffXStrehl(hw,  hi,  S)xT(hw,  hr,  S)  =  Id, 


where  T(hw,  hr,  S)  was  calculated  by  PcLnWin.  and  where 


I 


diff  - 


-f— f 

4UlJ 


P 


Specifically,  we  found  SI  and  S2  such  that 
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lditrxStrehl(hw,  hr,  S I  )xT(hw,  hi ,  S 1 )  <  Id 


and 


ldjfrxStrehl(hw,  hr,  S2)xT(hw,  hr,  S2)  >  Id 


By  continuity,  the  desired  value  of  S  lies  between  SI  and  S2,  and  we  determined 
it  by  bisection.  To  avoid  repeated  calculations  with  PcLnWin,  we  used  linear  in¬ 
terpolation  for  transmissivity  T,  after  SI  and  S2  became  sufficiently  close.  After 
bisection  converged,  we  compared  the  exact  value  of  T  for  the  converged  value  of 
S,  with  the  value  from  the  linear  interpolation  and  found,  typically,  3-  to  4-rigure 
agreement.  In  later  work,  we  used  the  analytic  models  of  PcLnWin's  e.xtinction  to 
determine  that  effect. 


Results 


The  leftmost  curve  of  Figure  14  shows  the  resulting  contour  of  I  kw/cm“  peak 
fluence  for  the  laser  described  in  [1 1].  That  laser  has  a  power  of  3  MW,  and  a 
mirror  diameter  of  1 .5  m. 

The  rightmost  contour  shows  the  diffraction  limit.  It  is  a  circle  of  radius  554  km, 
centered  on  the  weapon.  Please  note  that,  to  show  details  of  vertical  structure,  the 
chart  greatly  expands  the  scale  of  the  ordinate  axis  relative  to  that  of  the  abscissa. 
This  makes  the  circle  appear  straight  and  distorts  other  geometric  relations.  This 
contour’s  lowest  point  is  at  a  height  of  2  km,  the  lowest  point  for  which  the  re¬ 
fracted  ray  from  the  source  does  not  intersect  the  earth. 

Figure  14.  1  kW/cm^  Contours 


-  I  Surface  - Diffraction  Limit 

B  Weapon  — With  turbulence  and  extinction 
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Figure  1 5  shows  the  1  kw/cm“  contour  in  more  detail.  This  figure  suggests  that, 
even  without  benefit  of  adaptive  optics,  the  laser  described  in  the  APS  article 
could  deliver  damaging  fluence  to  a  target  despite  the  impacts  of  atmospheric  ex¬ 
tinction  and  turbulence,  at  tactically  useful  ranges. 

Nevertheless,  the  degree  to  which  such  a  laser  would  be  useful  as  a  tactical 
weapon  still  would  depend  on  effects  of  blooming  (which  we  consider  in  the  next 
subsection)  and  on  target  orientation,  geometry,  albedo,  and  motion  (which  we  do 
not  consider  in  this  report).  Jitter  in  the  airborne  laser  may  also  degrade  perform¬ 
ance;  we  consider  this  briefly  below. 

Figure  15.  Contour  Detail 


Distance,  km 


- Surface  ■  Weapon  -»—Wiith  Turbulence  and  Extinction 


For  insight  into  the  variation  of  the  range  of  a  laser  weapon  without  adaptive  op¬ 
tics,  Figure  1 6  shows  the  contour  of  0. 1  kw7cm^  intensity. 

Figure  16.  0.1  kW/cm^  Contour 


- Surface  B  Weapon  Wiith  Turbulence  and  Extinction 
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Blooming 


The  beam  heats  the  air  through  which  it  passes,  decreasing  refractive  index  un¬ 
evenly  and  distorting  the  beam.  This  effect,  called  “blooming,”  reduces  tluence  on 
the  target.  Here  we  consider  this  effect. 


Citing  an  article  by  Weichel  [8],  the  APS  report  says  that,  for  a  Gaussian  beam 
having  an  e-folding  radius  a(z)  where  z  is  path  length,  the  ratio  of  the  intensity  at 
the  target  with  blooming  to  that  without  blooming  is 


I(  Bloom)  I 

I(NoBIoom)  I-i-(BqIp)^ 


[Eq.  29] 


where 


diif, 

dT  4npoWoa^ 


[Eq.  30] 


and 


I  =  f 

'  L^Joa(z) 


aoWoa(z')To 

doa^(z')w(z>oT(z') 


J 


-dz' 


dz 


[Eq.3l] 


In  Equations  30  and  3 1 ,  a(z)  denotes  the  part  of  the  atmospheric  extinction  coefficient 
that  puts  power  into  the  air  through  which  the  beam  passes,  w(z)  the  effective  wind 
velocity  coming  from  the  motion  of  the  beam  relative  to  the  (possibly  moving)  atmos¬ 
phere,  T(z)  the  temperature,  and  a(z)  the  beam  radius,  all  varying  with  distance  z  along 
the  path.  The  subscript  0  on  any  of  these  denotes  the  corresponding  value  at  z  =  0. 


The  following  subsections  explain  how  we  evaluated  normalized  extinction  coef¬ 
ficient  beam  radius and  wind  profile  . 


ttr 


Wr 


EXTINCTION  AND  ABSORPTION  COEFFICIENTS 


Using  PcLnWin  outputs,  we  developed  analytic  models  of  the  variations  of  molecu¬ 
lar  absorption  and  atmospheric  extinction  coefficients  with  height. 

Figure  1 7  compares  our  model  of  molecular  absorption  with  the  one  used  hy  Walter 
and  Mani  [9]. 
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Figure  1 7.  Comparing  Molecular  Absorption 


- F^LnWIn  Data - WaKer  Model  A  LMI  Model 


Figure  1 8  compares  our  extinction  coefficient  with  the  one  used  by  Walter  and  Mani  [9], 
Figure  18.  Extinction  Coefficient  Model  and  Data 


- PcLnWin - Walter  a  Model 


Figures  1 7  and  1 8  show  that,  while  PcLnWin’s  molecular  absorption  agrees 
closely  with  Walter’s  results,  PcLnWin’s  total  extinction,  and  thus  by  implication 
aerosol  absorption  and  extinction,  are  considerably  larger  for  heights  of  2  to  5  km 
than  those  of  the  model  used  by  Walter  and  Mani. 

Beam  Radius 

For  the  variation  of  beam  radius  a(z),  we  assumed  linear  variation  from  the  mirror 
radius  at  z  =  0  to  the  radius  of  the  turbulence-broadened,  focused  spot  at  z  =  L.  By 
"turbulence-broadened,  focused  spot”  we  mean  what  Fante  [9]  calls  the  “short¬ 
term  beam  spread”  for  the  focused  beam.  As  Fante  explains,  long-scale  turbulence 
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causes  the  beam  to  wander,  while  short-scale  turbulence  spreads  it  out.  Thus  the 
target  receives  a  broadened,  wandering  spot  (Figure  19). 


Figure  19.  Broadened,  Wandering  Spot 

Y 


Source:  after  Fante  [9],  footnote  5. 

We  assume  that  the  heating  that  causes  blooming  is  influenced  by  the  broadening, 
but  not  by  the  wandering. 


Fante  gives  an  expression  for  the  radius  of  the  broadened  beam,  Ps 


4L^  4L- 

Ps  ~  ,  2t'«2  1  2  2 

^^0^  ^^Opf 


I  -  0.621  ^ 
D 


/3n6/5 


[Eq.  32] 


which  we  used.  The  parameter  pf  in  Equation  32  is  equal  to  48  percent  of  the  pa¬ 
rameter  ros  evaluated  in  Equation  22. 

Thus,  we  took 


a(z)  =  D-(D-p,).^. 


|Eq.33| 


which  leads  to 


^0  ^  D  J  L 


[Eq.  34] 


To  evaluate  the  iterated  integral  of  Equation  31  with  a  single  integral,  we  note  that 
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zioa(z) 


aoWoa(z')To 

■ioa-(z')w(z')aoT(z') 


J 


dz/ 


dz  =  j 


a5Woa(z')To 


=oa‘(z')w(z')aoT(z') 


-dz 


dz' 


Jl  a(z) 

[Eq.35] 


and  that  the  integral  with  respect  to  z  can  be  evaluated  in  closed  form,  as 


J^dz=  J _ ^ 

Aa{z)  JlI-(1-p,/D)z/L  I-p,/D 


In 


l-(l-p,/D)z'/L 

Ps/D 


=  g(z') 

[Eq. 36] 


With  the  g(z')  determined  in  Equation  36,  and  using  values  of  To/T(z')  from  the 
1976  U.S.  Standard  Atmosphere,  we  evaluated  Ib  numerically  as 


Ib  -  J 


aoWoa(z')To 
loa-(z')w(z')aoT(z') 


778(2')  dz' 


[Eq.  37] 


Wind  Profile 

We  needed  to  develop  a  wind  profile  w(z)  that  is  realistic,  and  for  which  the  inte¬ 
gral  of  Equation  37  converges,  and  for  which  the  Equation  D.29  of  the  APS  arti¬ 
cle  [  1 1  ]  will  apply.  That  equation,  specifically, 

pCpw|^  =  aI  [Eq.38] 

is  the  starting  point  of  the  blooming  analysis  of  that  reference. 

For  air-to-air  engagements  (except  for  shots  made  flying  directly  toward  the  target 
or  away  from  it)  the  wind  relative  to  the  beam  will  be  large  along  the  entire  path. 
Equation  38  will  apply,  and  the  integral  will  certainly  converge. 

But  for  a  stationary  target  on  the  ground,  it  will  not  do  simply  to  take  the  velocity 
profile  generated  by  the  motion  of  the  beam  from  aircraft  to  a  fixed  ground  target. 
That  wind  profile  may  fall  linearly  to  zero  at  the  target,  and  the  integral  of  Equa¬ 
tion  37  would  not  converge. 

It  is,  however,  possible  that  typical  winds  would  be  large  enough  for  Equation  38 
to  apply.  If  so,  then  a  simple  profile  w(z)  that  began  with  a  speed  like  that  of  the 
weapon-bearing  aircraft  at  z  =  0,  and  fell  to  a  typical  surface  wind  speed  at  z  =  L, 
could  be  used. 

We  began  our  consideration  of  wind  profile  with  the  equation 
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[Eq.  39] 


— +  q  VT  =  — V-T  +  -!-Q. 

pCp  pCp 

expressing  conservation  of  energy  in  an  isobaric  motion  of  a  tluid  with  constant 
specific  heat  at  constant  pressure  Cp  and  thermal  conductivity  k.  neglecting  vis¬ 
cous  dissipation.  Equation  39  is  written  for  an  inertial  frame  moving  with  the 
beam  locally.  In  Equation  39, q  is  fluid  velocity  in  cm/sec,  T  is  temperature  in  “K, 

Q  denotes  rate  of  heat  addition  in  watts/cm^,  k  is  thermal  conductivity  in 
watts/(cm  °K),  and  p  is  density  in  g/cm^. 

Introducing  non-dimensional  variables  T,  q,  x  and  twith  the  definitions 


pcpd^ 
q  =  q/w 

X  =  X  /  d 

t  =  T/T,,f 


[Eq.  40] 


where  x  is  position  coordinate,  w  is  tluid  velocity  relative  to  the  beam,  and  d  is 
beam  diameter,  brings  Equation  39  to 


wpc  d  3t 


9T  _  „ 

+  q . VT  = 


wpcpd 


V‘T  + 


w 


dT, 


Q 


ref 


[Eq.4l] 


In  Equation  41  we  have  suppressed  the  ^  symbols  for  clarity.  Equation  40  shows 
that  temperature  changes  will  satisfy 


q  VT  = 


lEq.42] 


provided  that - «  I ,  or,  equivalently,  that  w  » - .  For  a  representa- 

wpCpd  pcpd 

tive  case,  k  =  2.55x  1 0“*  joule/(sec  cm  "K),  p  =  1 .225x  1 0'^  g/cm^,  Cp  =  I  joule/(g 
°K),  and  the  condition  that  Equation  41  governs  temperature  changes,  and  is  w  » 
0.01  cm/sec. 

This  condition  is  likely  to  be  met  by  ordinary  wind,  even  if  the  beam  is  fixed  in 
space,  as  it  would  be  on  a  stationary  target  on  the  ground. 

With  this  result  in  mind,  we  developed  the  profile  of  wind  relative  to  the  beam  in 
this  way:  If  Vpand  Vj  denote  respectively  the  velocity  vectors  of  platform  and 
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target,  and  if  W  denotes  the  wind  velocity,  then  at  a  point  on  the  beam  located  a 
fraction  q  of  the  distance  from  platform  to  target,  the  beam  will  have  velocity 

Pig)  =  V,  +  giV,  -V,)  +  IV ig))  [Eq.  43] 


The  magnitude  Vpcrp  of  the  velocity  perpendicular  to  the  beam  will  be  given  by 
the  magnitude  of  the  cross  product  of  a  unit  vector  parallel  to  the  beam,  and  the 

beam's  velocity  V .  We  evaluated  that  cross  product,  arbitrarily  choosing  the 

components  of  W  to  be  a  constant  2  m/s  in  each  of  3  perpendicular  directions,  to 
obtain  w(z). 


In  addition  to  parameters  determined  in  evaluating  1b,  evaluating  Bo  required  only 

a  value  for  —  .  We  used  the  value  -Ixl  0'^  cmV(g  ”K)  given  in  the  APS  re- 

Po  dT 

port  cited  above. 


Blooming  Results 


With  the  approach  described  above,  we  found  that  blooming  would  reduce  the 
range  at  which  a  surface  target  receives  I  kW  per  square  centimeter  from  71  km 
to  about  58  km,  a  reduction  of  roughly  18  percent.  In  the  case  of  0.1  kw/cm“, 
blooming  reduced  surface  range  only  slightly,  from  1 13  km  to  1 12  km.  The  rea¬ 
son  for  this  is  that  turbulence  causes  a  larger  spot  size  in  the  0.1  kW/cm“  case  than 
it  did  in  the  case  of  1  kW/cm2.  There  is  less  heating  in  the  less-concentrated 
beam.  and.  consequently,  less  blooming. 

Stable  Secondary  Atmospheric  Motions 

Long-range  propagation  through  the  atmosphere  is  certain  to  encounter  turbu¬ 
lence.  Also,  as  even  casual  observation  of  the  sky  will  verify,  stable,  wave-like 
motions  of  the  atmosphere  are  commonplace  (Figure  20). 
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Figure  20.  Atmospheric  Waves,  Made  Visible  by  Clouds 


We  made  only  a  preliminary  consideration  of  the  potential  effects  of  atmospheric 
waves  on  airborne  laser  weapons.  If  (but  only  if)  it  is  accurate  to  regard  the  waves 
as  ‘‘frozen.”  that  is,  varying  in  time  and  space  so  slowly  that  our  refracted-ray  cal¬ 
culations  give  useful  indications  of  the  waves’  effects  on  the  laser  beam,  then 
those  calculations  may  be  used  to  study  the  displacement  of  the  beam  by  atmos¬ 
pheric  waves.  We  carried  out  a  few  such  calculations,  using  the  atmospheric  wave 
model  given  by  Wuertele  et  alii  [  12].  The  results  indicate  that,  for  propagations 
over  ranges  of  the  order  of  100  km,  atmospheric  waves  might  result  in  beam  dis¬ 
placements  on  the  order  of  one  centimeter. 

We  believe  that  the  effects  of  atmospheric  waves  on  airborne  laser  weapons  de¬ 
serves  more  study  than  we  were  able  to  give  it  in  this  project. 
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Summary  and  Conclusions 


Significant  progress  has  been  made  during  this  year  in  both  combat  modeling  and 
laser  propagation  modeling.  The  air-to-air  and  air-to-ground  campaign  analysis 
has  addressed  several  ad  hoc  problems  and  has  reached  the  point  where  the  cam¬ 
paign  model  can  be  used  for  TACAIR  analyses  and  decision  making. 

The  air  superiority  campaign  model,  SLAACM,  has  been  incorporated  into  an  MS 
Excel  /  Visual  Basic  program.  The  current  version  of  SLAACM  models  4-ship  Blue 
defense  formations  versus  12-ship  Red  attack  packages.  The  Red  packages  include 
two  4-ship  formations  of  escorts,  protecting  a  4-ship  formation  of  bombers.  There 
are  quality  options  for  both  Blue  and  Red  fighters  and  for  Red  bombers. 

For  each  day  of  the  campaign.  Red  assembles  its  attack  packages  to  maximize  its 
bombing  effectiveness  and  the  destruction  of  Blue  forces.  As  a  model  option. 

Blue  can  respond  randomly  to  the  Red  attack  or  use  prior  knowledge  to  preferen¬ 
tially  attack  the  highest  threat  Red  packages. 

The  model  output  includes 

♦  updated  daily  listings  of  the  Red  order  of  battle, 

♦  graphical  depiction  of  Blue  and  Red  fighter  and  Red  bomber  inventories,  and 

♦  tabulated  values  of  targets  destroyed  and  the  tons  of  bombs  delivered. 

A  campaign  of  several  days  runs  in  less  than  I  minute  on  a  Pentium  3  PC. 

The  laser  analysis  has  reached  the  point  where  the  propagation  modeling  is  well  un¬ 
derstood  and  efficiently  modeled.  Addition  of  target  interaction  and  scaling  informa¬ 
tion  related  to  high-energy  laser  systems  must  be  addressed  in  future  work  to  enable 
comparison  of  high-energy  lasers  with  current  air-to-air  and  air-to-ground  weapons. 

The  laser  propagation  algorithms  are  incorporated  in  an  MS  Excel,  Visual  Basic  work¬ 
book  adequate  for  analysis,  but  have  not  yet  been  structured  into  a  publishable  model. 

Our  current  version  models  the  1 .3 1 5  |i  iodine  laser  with  a  limited  number  of  atmos¬ 
pheric  environments.  The  model  runs  quickly  and  produces  tabular  laser  intensity  and 
geometry  data.  Graphical  data  are  also  produced  but  are  not  updated  run  by  run. 

We  plan  to  expand  tbe  model  to  include  a  solid  state  laser  wavelength  and  additional  at¬ 
mospheric  molecular  and  aerosol  environments.  Following  clean-up  of  coding,  implemen¬ 
tation  of  improved  input/output  features  and  additional  verification  and  validation,  we 
intend  to  publish  the  model  as  the  Airborne  Laser  Infrared  Transmission  (ALIRT)  model. 

Finally,  we  plan  to  investigate  the  use  of  the  methods  developed  in  the  laser 
propagation  analysis  to  investigate  the  performance  requirements  for  laser  range¬ 
finders,  laser  designators,  and  electro-optical  sensors  both  at  independent  sensors, 
and  as  part  of  a  high-energy  laser  weapon  system. 
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